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ABSTRACT 


l  > 

'packaging^ 


This  report  shovs  systematic  techniques  for  handling  or  "''packaging' 
of  n- tuples  of  n- tuples  of  position,  velocity^and  acceleration  vectors . 
large  systems  of  vectors  cast  Into  stete-vector  representations  lead^  In 
a  natural  way^  to  natrlx  Riecatl ^aquations. 
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The  scalar  Rlccati  Equation  as  given  by  Ince  (Reference  l)  is 
y  +  ay8  +  by  +  c  »  0  (i-l) 

where  y  is  a  scalar- valued  variable  of  the  scalar  argument  (variable) 
tine,  a,  b,  and  c  are  constants  pr  tine  varying  coefficients. 

The  matrix  Rlccati  Equation  as  given  by  Levine  (Reference  2)  is 

* 

* 

t+BiY  +  IB#  +  YAY+C-(0]  (i-2) 

where  Y  (t)  is  a  matrix- valued  function  of  a  scalar  argument  (single  variable 
tiae)  and  A,  Bj,,  Bg  and  C  are  constant  or  time-varying  coefficient-matrices  * 
Kalman,  Bellman,  Fontrjagin  and  other  researchers  working  in  the  area 
of  mathematical  programing  and  optimisation  frequently  refer  to  the  matrix 
Rlccati  Equation.  The  detemlnatlon  of  an  optimal  control  In  feedback  control 
theory  or  of  an  optimal  filter  in  estimation  and  prediction  theory  invokes 
a  consideration  of  the  Rlccati  Equation. 

This  report  shows  how  numerously  one  can  obtain  matrix  Rlccati  relations 
in  vector  dynamics.  A  geometrical  insight  into  the  relations  described  by  the 
matrices  is  presented. 
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Section  I.  Some  Aspects  of  Generalised  Qrnamlcal  Systems 

This  section  considers  systems  of  noters  and  their  first  and  second 
derivatives.  Representations  of  the  vectors  independent  of  their  coordinates 
in  bases  are  presented  parallel  to  matrices  cf  their  co ordinate b  relative  to 
movirg  bases.  Metrics  of  invariant  inner-products  of  velocity  vectors, 
representing  a  matrix  of  scalars  of  the  motion  of  a  system  of  vectors, 
motivates  the  notion  of  generalised  energies.  A  scalar  invariant  of  the 
system  is  generated  as  the  trace  of  the  matrix  and  has  lte  analog  in  total 
system  energy. 

Many  other  classical  mechanics  notions  are  presented  in  a  unifying 
systematic  manner  based  on  fundamentals  of  mathematical  transformation  theory. 

Consider  a  vector  x  in  real  n- space  with  coordinates  in  a  fixed  basis 


<f  and  in  two  other  arbitrary  bases  <^y  and <^z,  that  is 

x  ■  x^  .  (1) 

The  coordinates  of  the  <^y  basis  in  the^f  space  is 

<y  -<?  *f  (2) 

and  similiarly 

<7-<*Z*  (3) 

The  system  of  2n  +  1  vectors 

<^s  -  (x,<y,<I)  »  (x,  yi,  y2,  y3,  *i>  za)  00 


in  n-  space  with  no  constraints  represents  n(2n  +  1)  independent  coordinates 
that  is  2na  +  n  scalars  (field  elements )  specify  the  configuration  of  the 

90 

seven  vectors.  Kach  second-derivative  vector  (e.g.  x  in  n- apace  requires 
2  n  Integrators  (vlth  their  initial  conditions)  to  obtain  the  vector. 
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Hence  the  dynamical  propagation  of  the  system  of  vectors 

"  (*#^y>^*)  -  (x,  yi,  y8,  ya,  t\,  za,  z3)  (5) 

requires  2  x  n  x  7  *  14  n  integrators.  In  three  space  one  has  42  integrators 
in  a  computer  mechanization  of  the  system.  The  time  derivatives  of  (?)  are 
with  respect  to  the  fixed  background  basis,  that  is 

<f -<3  •  (6) 

Using  (2)  and  (3)  in  (l) ,  one  obtains 

x^>  ■  Yf  x*>  »  Z*  x^>  (7) 

hence  for  oomplete  flexibility  in  mapping  between  bases,  the  transformation 
matrices  I  and  Z,  must  be  obtained.  These  transformations  as  well  as  a 
number  of  other  relations  are  obtained  as  solutions  of  Matrix  Riceatl  type 
equations. 

Before  bogging  down  in  the  detailed  maze  of  matrices  of  coordinates, 
let  us  glimpse  at  the  overall  second-derivative  picture  of  the  system  of 
vectors.  Bote  that  above  we  state  matrices  of  coordinates.  This  handling 
of  large  packages  of  coordinates  (matrices)  alleviates  the  tediousness  of 
the  standard  classical  treatise  of  mechanics  systems  of  particles  wherein 
one  deals  with  a  maze  of  super  and  subscripts  which  rapidly  swamp  the 
meaning  of  the  physical  picture. 

A.  Matrices  of  Invariant  Representations 

The  relative  configuration  of  the  system  of  vectors  may  be 
described  by  a  matrix  of  inner-products  (independent  of  coordinates)  by 
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t>  < 


8  m 


or 


Qss 


X.X 

>•<? 

x* 

y>.x 

y><y 

£>•<« 

^  .X 

*><*> 

where  geometrically  one  has  the  conventional  notion  of 
X  .  X  -  |x|  |x| 

x  «  (yx*x,  ya»x,  ya«x) 

y^>  •  < 


■t 

yi**i 

yi**2 

yx**a 

ya-*i 

ya**a 

ya-*a 

ys**i 

ys*« 

ys**s 

Sr*  » 


One  nay  now  write  (9)  as 


Qs8  j 

*•<7 

x.<^i 

5>-5 

Sry 

Sr* 

Qtz 

The  first  derivative  dynamic*  of  the  ays  tarn  la  daacribad  hy 

0,8  -  B>  <1  +?><S 


(8) 


(9) 


(10) 

(11) 

(12) 


(13) 


(14) 


or 
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Qbs 


/ *  \  (*»</#  <C*)  V*  V  (x,<Cy»<$) 


(15) 


y> 


Bote  that  the  tine  derivative  terms  of  (15)  are  coordinate  free.  In 
later  sections  the  matrices  of  coordinates  viU  enter  the  picture. 

The  second  derivative  dynamics  of  the  system  is  described  by 

Qss  -^>*<^8  +  2  s>*<8  +  (16) 

Clearly  one  detects  the  notion  of  total  system  Kinetic  energy  or 
rather  a  matrix  of  some  of  the  scalars  of  the  motion  in  the  matrix  of 
inner-products  i>  <i,  that  is 


i<7 

t><5 

t>  <3 

The  total  "energy"  is  the  trace  of/ (17),  that  is 
e  «<^8*8^>  »  trace  ^ 


(17) 


(18) 


or 


<ii>  ■  x.x  +  trace  <Cy  +  trace  i><i 

—  -  Zk  >*  'v 

e  -  x«x  +<^y*y> +<J*  V>  •  (19) 

The  trace  of  a  matrix  is  an  invariant  (independent  of  bases. ) 

The  analogs  of  the  classical  concepts  of  energy,  power,  etc, lor  a 

system  of  vectors  may  also  be  formulated  and  handled  at  the  matrix  of 

vectors  level.  Lator  reports  will  present  these  generalized  systems  of 

vectors  in  generalized  coordinates  in  a  Hamiltonian  framework. 
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B.  Matrices  of  Matrices  of  Coordinates 

Consider  the  matrix  of  coordinates  of  the  ^  vectors  in  the^^y 
basis,  that  is  by  (2)  and  (3) 


<cf  -  Y-Iz  -  T  . 
The  system  of  vectors  is 

\s  -  (x,<^y,<^y  t) 


(20) 


now 


(21) 


and 


t><* 


x,x 


X*  </  X*  <F* 

V  *VyT 


^>-X 

T'£>  -  T\y  T 


(22) 


The  matrix  of  scalars  T  of  (20)  may  be  evaluated  as  a  matrix  of  inner- 
products  by  projecting  the  vectore  onto  the  dual  vectors  ^  ,  that  is 

f><5  -  [?i  *j]  - 1^"<?  *  -  *  > 


(23) 


since  by  definition  of  duals  or  reciprocals 


0><y 


I  . 


(24) 


It  is  shown  in  later  sections  that  the  matrix  T  satisfies  a  matrix 
Rlccati  differential  equation. 

C.  The  Allas  -  Alibi  Aspect 

In  algebraic  transformation  theory  one  is  confronted  with  the  vector 
matrix  equation 

w^>  -  A^>  (25) 
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which  has  two  interpretations; 

(1)  and  are  two  different  vectors. 

(2)  and  are  the  sane  vectors  in  different  coordinate 

systems  (bases).  Clearly  (25)  is  the  same  form  as  (7). 

Halmos*points  out  that  the  theory  of  changing  bases  is  coextensive 
with  the  theory  of  Invertible  linear  transformations .  Ha  states  that  an 
invertible  linear  transformation  is  an  automorphism,  whereby  an  auto¬ 
morphism  means  an  isomorphism  of  a  vector  space  with  itself.  Balnos 
points  out  that  conversely,  every  autostorphlsm  is  an  invertible  linear 
transformation . 

The  study  of  the  properties  of  a  linear  transformation  A  on  a  vector 
“>  and  of  the  structure  of  the  transformation  A  necessitates  that  one  study 
the  Effect  of  A  on  n- linearly  independent  vectors  i, 

2*  •  ••  u>n- 

Jtor  example,  in  the  three  space  three  linearly  independent  vectors 
u>  ±  in  the  domain  space  under  A  map  to  three  vectors  in  the  Range 
space  with  respect  to  A  as 

^  3  ]  “  [A  A  A  a]  ^ 

or 

W  ■  A  U  (27) 

Since  U  are  linearly  Independent  their  matrix  of  coordinates  U  is  an 
invertible  matrix,  hence 

-1 

A  -  W  U  .  (28) 


*Halmos,  P.  86. 
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Clearly,  then  in  the  study  of  the  dynamios  of  the  single  rector  x, 
one  must  study  the  faction  of  n  linearly  independent  rectors 

-  fa. $  *8  ...  xn)  (29) 

In  differential  equation  theory  this  is  equivalent  to  obtaining  n 
linearly  independent  solutions. 

For  example,  if  the  vector  x  is  known  to  propagate  under  a  linear 
lav  Aft)  where  A  is  a 

x>  -  A(t)x>  (30) 

function  of  time,  one  must  consider  the  extended  system  of  (4)  as 

(31) 

to  obtain  tho  dynamics  of  x  In  tvo  arbitrary  cats  of  moving  bass.  «  and 
<?  In  other  words,  one  must  consider  n- linearly  Independent  solution 
vectors. 

This  paper  will  not  go  into  the  study  of  the  analytical  properties 
of  the  systems  of  dynamical  vectors  of  (31). 
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Section  II.  Vector  System  Equations  of  l^btion 

This  report  is  Halted  to  formulation  of  the  equation  of  notion  for 
the  system  of  (l«4)  for  computer  mechanization  and  to  yelld  geometrical 
insight  into  the  geometrical  dynamics. 

A.  First  Derivative  Vector  Matrix  Relations 

The  rector  x  by  (£-1)  in  the  two  arbitrary  spring  bases  is 

x  -  ^  ^  x*>  -  (i) 

and  the  derlratlre  with  respect  to  time  is 

5-<^z3^>+<3r»h>  -<5x^>  +<tx5>.  (2) 

Since  <y  is  a  basis,  the  three  Telocity  vectors may  be  expressed 
in  the  <^y  basis,  that  is  a  matrix  of  scalars  exist  such  that 


<5  -<y  v  <y) 

and  similarly  for  the  rectors 


(5) 


By  (1-2) 

<5  *<^i  V  (z) 

and 

<?  "<?* 

^  -1 

(5) 

hence 

<7  -<?* 

(6) 

<y  .<?  I 

(7) 
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since  by  (1-6) 


<f.<(5  (8) 

hence 

.j,  "l* 

<y  -<y  y  t  .  (9) 

By  (3)  and  (9) 

-1. 

Y  Y  -  V  (y)  ,  (10) 

and  similarly  for  the  basis  ^z,  that  is 

.  -1  . 

<1  -<Z  Z  Z  -<S  V  (z)  (11) 

By  (10)  and  (ll) 

(12) 
(13) 

The  two  first  order  matrix  differential  equations  above  are  special 
cases  of  the  Rlccatl  matrix  equation  shown  ln(l-2). 

By  (1-20) 

<*  t  m<y  y  1  z  d4) 

hence  also 

<*  -<c y  t  +<o  *  (15) 
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Using  (9)  and  (ll)  in  (15) 


T  ~<[y  T  V  (2)  +<^y  V  (y)  T 


(16) 


or 


T  -  T  V  (2)  +  V  (y)  T 


(17) 


which  is  also  a  matrix  Riccati  equation  as  seen  by  (1-2). 

One  may  derive  (17)  directly  from  the  matrix  relations  of  (14),  that  is 


-1 

T  -  Y  Z 


(18) 


and 


now 


hence 


therefore 


# 

T  - 


-1  -1  . 

y  ' z+Y  z 


-1 

Y  Y  -  I 


-1  -1  . 

d  Y  (Y)  +  Y  Y  -  [0 ] 

dt 


-1  -1  .  -1 

d  (Y  )  -  Y  Y  Y 

dt 


(19) 


(20) 


(21) 


(22) 


Equation  (22)  in  (19)  yields 


.  -1  .  -1  -1  . 

T-Y  YY  Z+Y  Z 


(23) 
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Using  the  identity 

-1 

Z  Z  -  I  ,  (24) 


(23)  becomes 

•  -1  •  "1  “1  -1  • 

T  -  Y  Y  (Y  Z)  +  (Y  Z)  Z  Z  (25) 


Using  (18)  in  (25) 

-1  .  -1  . 

T-Y  YT  +  TZ  Z  (26) 


Clearly  if  Y  and  Z  satisfy 

Y  -  Y  V  (y)  (27) 


and 


Z  -  Z  V  (*)  (28) 


then  (26)  may  be  written  as 

T  ■  V  (y)  T  +  I  V  (1)  (29) 

which  agrees  with  (17)* 

The  latter  derivation  is  mechanical  manipulation  and  does  not  convey 
the  direct  geometrical  clarity  as  does  the  first  derivation  leading  to  (17)* 
Consider  the  metric- natrix  of  the  ^  basis  N yy,  that  Is 

*Vy  -  y>*  <y  -  I'Mff  Y  <30) 
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The  derivative  of  (30)  is 

%  -  O*'  <^y  +  y^>  <Cj  •  (51) 

By  (3) 

<y  ■  <^y  V  (y)  (32) 

and  transposing 

*  V  (y)  y^>  (33) 

hence  using  (32)  and  (33)  in  (31 )  yields 

(3*0 

which  once  again  is  a  special  case  of  the  Matrix  RLccatl  equation  (i-2). 

A  repetition  of  the  same  ar gurnet  for  the  dynamics  of  the  z 
basis  yields 

(35) 

Consider  next  the  matrix  of  inner-products  between  the  bases 
<^y  and<^z,  that  is 

Q  -  y^>  <^z  -  Hyy  T  (36) 

The  matrix  Q  specifies  the  relative  orientations  of  the  vectors. 

The  derivative  of  (36)  is 

i-y><* +y>*<*  (37) 
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Using  (33)  and  (ll)  in  (37) 

Q  -  V \y)  y^>’<ji  +  (*)  (38) 

and  by  (36) 


B.  Second  Derivative  Vector-Matrix  Relations 

The  vector  x  in  the^?,^y  and  the<^s  bases  by  (l-l)  is 


x  -<^J?  «^y  x*]>  (40) 

and  the  velocity  vector  is 

5  •<?  i£>  -<y  <M) 

»^s  x*>  +  ^  x^  .  (i»2) 

The  acceleration  vector  is 

x  •<7x1^«<j^>  +  2^/>+<j5r^>  (43) 

»<^t  x^  +  2  x^  +^5  x^  .  (W) 
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Section  III.  Trajectory  of  a  Vector  x  and  one  Arbitrary  Time  Varying  Basis  <^y. 

The  first  case  vill  consider  the  simple  system  of  four  arbitrary 
vectors  with  the  only  condition  that  the  three  vectors  ^y  fofta  a  basis, 
that  is 

<J  -  (x,<^y)  -  (x  (t),  yi(t),  y2(t),  y3(t))  (45) 

vhere 

x  -  <?  xf(t)>  -<y  xfr>  (W) 

and 

<7  »<^?  Y  (t)  .  (47) 

In  three  space  we  have  a  system  of  four  vectors  each  having  three 
independent  variables  or  a  total  of  twelve  arbitrary  coordinates  to 
specify  the  geometrical  configuration  of  the  system  of  vectors. 

Clearly,  at  the  velocity  level  we  need  twelve  integrations  and 

<*  -<?  [,£>  -  iw]  -<?  m 

5x4 

twelve  initial  conditions  or  constants  of  the  motion  as  shown  by  (48)  where 
the  components  of  the  3x4  matrix  are  all  in  the  "fixed* bases  ^^7. 

In  most  real-life  problems  the  parameters  of  the  system  and  the 
constants  of  the  system  occur  most  naturally  in  moving  bases. 
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Clearly,  the  second  derivative  vectors  of  the  system  require  24 
Integrations,  that  is 

<s  -<^7  ,  Y(t)  J  -<v  (49) 

and  24  initdl  conditions,  the  additional  twelve  are  conditions  oz> 
Initial  velocities,  since 


J >  *  v(x)£> 

and 

Y  -  V(y) 

hence 

<T®  -<7  [4(x)£>  ,  V(y)  J  (50) 

♦t  ^ 

The  acceleration  vector  x  In  the  moving  basis  ^y  by  (45)  is 

x  x^  +  2^y  +^y  x1^  .  (51) 


By  (32) 

<5  -<y  V(y) 

and 

"Cy  m<Cy  V(y)  +<^y  Y(y) 
or 

<Cy  [v*(y)  ♦  V(y)  J 


(52) 


(53) 

(54) 


16 


Using  (52)  and  (55)  in  (50) 

x  «^y  +  2  V(y)  x^  +  £v*(y)  +  V(y)  j  xJ^  J  .  (55) 

also 

x-3.-<^ya^  (56) 

vheregj^msy  be  considared  the  driving  vector  (specific  force  in  Mechanics), 
hence 

x£>  +  2  V(y)  x3^>  +  (V*  +  V)  x3^  -  a/>  (57) 

Equation  (57)  is  a  second-order  matrix  differential  equation  linear  in 
the  vector  variable  xJ^>  and  its  derivatives,  that  is 

V>  +  Ae(t)  3^  +  Ai(t)  x>  -  a^  (58) 

vhere  in  general  the  vector  ^  is  a  non-linear  function  of  x]>  and 
*>  ,  that  is 

a>  -  a  (x>  ,  (59) 

depending  on  the  nature  of  . 

It  has  been  observed  that  we  need  24  integrations  for  the  systems 
under  consideration*  Six  integrations  are  required  for  (58),  however, 
the  position  and  the  velocity  coefficient  matrices  Ai  and  Ag  must  be 
determined.  Observe  that  one  vector  x  requires  six  Integrators  and  that 
the  system  of  three  vectors  <^y  will  require  18  more  integrations. 
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By  equation  (5*0 


<0  m<Cj  [*®(y)  +  v(y)j  "<3r  H 


(60) 


that  is 


(yi *  Ya>  y3)  -  (hi,  Ua,  ha)  (6l) 

and  the  three  "driving"  vectors  <^5  have  coordinates  in  the  <^y  basis  as 

<J  -<y  H  .  (62) 

Thus,  if  the  matrix  H  is  knovn  the  first  order  non-linear  matrix 
equation  can  be  mechanized,  that  is 


V(y)  +  V*(y)  -  H  . 


(63) 


Equation  (63)  requires  nine  Integrations  and  nine  initial  conditions. 
Equation  (27)  provides  the  other  nine  Integrations,  that  is 


T  -  I  V(y) 

The  system  of  vector  and  matrix  differential  equations  is 


(64) 


(65) 
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If  the  constraints  are  Imposed  on  the  basis  <^y,  that  they  are 
orthonormal,  that  is,  that 

Myy  -  y>  <y  -  1  -  [yi-yj]  (66) 

then  the  symmetric  matrix  contains  six  independent  conditions  to  be 
satisfied.  These  conditions  are 

yi*yi  -  1  i  ■  1,2,3  (6?) 

vhich  is  three  conditions,  and 

yi-yj  -  0  i  \  j  (68) 

since 

yi-yj  -  yyyi  (69) 

the  six  off-diagonal  conditions  are  redundant  because  of  symmetry,  that 
is  only  three  independent  conditions. 

From  the  above  discussion,  it  is  clear  that  the  three  rectors  <^y 
with  the  constraint  (66)  are  completely  specified  by  three  independent 
coordinates  instead  of  nine,  that  is  the  nine  elements  of  the  matrix  Y 
are  functions  of  only  three  variables  q^>  ,  hence 

<^y  Y  (<U,  92,  93)  •  (70) 

It  is  well  known  that  when  <^t  is  also  op  O.V  basis,  or 

«ff  -I 

p>  -  Y'f> 
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(71) 

(72) 


hence 


y>*<y-i'MffY  .  (73) 

By  (66)  and  (7l) 

I  -  t'i  (74) 

hence 

/ 

Y  «  Y 

and  T  is  said  to  be  an  orthogonal  matrix. 

The  theory  of  Suler  angles  and  relations  to  0.1  bases  is  veil  known,  that  is 

<y  -<*  *(♦»)  «)(«*)  *i(*i)  -<**(♦>  ) 

where  are  luler  angles  and  are  rotation  matrices. 

It  will  now  be  shown  that  the  matrix  Y  (y)  of  (31)  is  skew- symmetric. 
Consider 

y>>  *  <5  -  i 


and  differentiating 

t>  •  <y  +  £>,<0  -  l°) 


hence 


y><y 


t><5 


(75) 
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Transposing  (51) 


•  V  '(y)  y^> 

projecting  (76)  ont-'  <^y,  that  is 

«  nl  (y)  r>  •  <^y  -  V^y) 

and  projecting  (51)  onto  y^>  that  \r 

y>>  •  <C?  -  v(y) 

using  (77)  and  (78)  in  (75) 

V  '(y)  -  -  V(y) 
or 

v(y)  -  -  ?'( y) 


(76) 


(77) 


(78) 


(79) 


hence  the  matrix  V  is  sk.ev-  symmetric .  Define  its  el— nts  as 


V 


-life  0 


(80) 


The  relationship  to  Buler  angles  is  given  in  reference 
The  matrix  equation  (65)  becomes 


Equating  the  indicated  elements  of  Equation  (8l) 


Wj.  +  (ife  ■  has 

wa  +  ^  «i  -  h*!  (82) 

Ufe  +  Wi  Wg  ■  hl2 

Equation  (82)  may  be  put  into  matrix  font  as 


The  above  non-linear  transformation  matrix  is  singular  (non- invertible) 
and  is. not  unique,  for  example 


Observe  that  the  transformation  matrix  of  Equation  (85)  is  a  function 
of  k>i  and  u<3  vhereas  the  matrix  of  (84)  is  a  function  of  «i  and  Ug. 

If  we  make  (83)  look  like  a  vector  in  the  <5  space 


<y  w>  +  <y  t  (w> )  w>  «<^y 
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then  the  vector  components  are 


h»a  *  ya'ks 

hgi  ■  ya*^i 
hi2  -yx*:52 

and 


Clearly  (63)  requires  only  three  integrations  rather  than  nine  as  before. 
Hence  the  system  (63)  reduces  to  18  integrators.  However,  these  are  not 
linearly  Independent  differential  equations  as  is  knpvn,  only  twelve  are 
llneaxly  independent.  Let  us  find  the  redundant  ones. 

Consider  Equation  (64),  that  is 

V  ($> )  -  Y'1  i 

Equating  the  three  unique  elements  of  V  in  (63)  to  the  right  hand 
matrix  product  yields 

wi  ■  wi  ^  ^ 

■  ufe  ^  (05) 
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The  relations  (85)  are  veil  known  to  imply 


w>  -T  (>)  ♦> 

where  T  )  Is  invertible  so  long  as  the  second  angle  is  not  90°,  that  is 

$>  -  T'1  «>  .  (86) 

One  can  nov  use  three  integrations  to  obtain  ^  and  then  generate 
the  trigonometric  relations  of  the  matrix  Y  since 


Y  -  I  ($» 


The  constrained  system  nov  has  twelve  integrations 

x 7^  +  2  V(t^)  +  ^Y(w)  +  V*)  x^  ■ 

+  T(^  )  ■  h£> 

£>  -  T‘1(#)  w]> 


(87) 


The  above  system  of  equations  looks  like  the  standard  six  degree  of 
freedom  rigid-body  dynamics.  Clearly  twelve  integrations  axe  required. 

The  transfer-matrix  flow  diagram  of  the  system  (87)  is  shown  in 
Figure  (l). 
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Section  IV  ,  The  Trajectory  of  a  Vector  x  and  tyo  Arbitrary  Time-Varying 
Bases <^y  and^l  (42  Integration*). 

In  aerospace  vehicle  math-modelling  one  must  consider  many  time 
varying  bases.  This  section  considers  the  system  of  vectors 

(x,  <Cy,  <$)  M 


and  the  relations 


x  x1^  -^r 

(2) 

vhere 

*3 

• 

V 

(5) 

or  by 

<7  -<7  i'1  z  -<y  T  . 

(4) 

Z  -  T  T 

(5) 

By  analogy  with  (II-60) 

[i-w  +  v(*)J 

(6) 
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The  system  eq.uat.ions  are  by  (11-65)  (11-17)  5 


*x£>  +  2  V(y)  +  £*(y)  +  V(y)J  ■  a^ 

6  integratora 

• 

V(y)  +  V*(y)  ■  H(y) 

9  Integrators 

V(z)  +  V*(z)  -  H(z) 

9  integrators 

T  +  Y(y)  T  -  T  V(z)  -  [0] 

9  Integrators 

Y  ■  Y  T(y) 

Algebraic  Relations 

Z  -  Y  T 
x*  >  -  Y  x*> 

and 

x*>  -  T  1  x7> 

9  Integrators 

The  above  system  requires  42  integrations  plus  the  indicated  algebraic 
manipulations.  The  system  is  not  unique,  that  is  one  could  obtain  many  other 
math  models  of  the  same  dynamical  system.  Jbr  example,  oonslder 

»<3r  T 
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and 


<5  * +<J  * 

and 

<C*  «<^y  T  +  2<J  T  +<5  T  »<5  H(*) 
by  (II-3)  and  (II-60) 

<C*  -<(/  jv(y)  +  2  Y“(y)j  *  +  T  +<C/  ? 


T  +  2  V(y)  T  +  |^V(y)  +  2  7*(y)J  T  -  H(t)  .  (8) 

Hence  the  system  is  defined  by 


the  transformation  matrix  I  is  a  function  of  only  thrae  parameters  -  —riy 
the  Buler  angles  orienting  the  two  bases  then  one  redness  tbs  system  (9) 
to  a  system  of  9  degrees  of  freedom  or  l8  integrations*  The  details  of 
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the  derivation  are  not  presented  In  this  paper.  The  classical  approach 
to  this  case  Introduces  such  "poor"  notions  as  "pseudo  vector"  or  axial 
vectors  and  relations  like 

*  -  * 7 

where  the  "axial"  angular  velocity  vectors  and  ay  are  "Inertial  angular 
velocities"  and  Ugy  is  the  relative  angular  velocity  vector. 
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